WELL-POSEDNESS OF NEMATIC LIQUID CRYSTAL FLOW IN L 3 j 

JAY LAWRENCE HINEMAN AND CHANGYOU WANG 

Abstract. In this paper, we establish the local well-posedness for the Cauchy problem of the simplified 
version of hydrodynamic flow of nematic liquid crystals ll.lt in M 3 for any initial data («o> <^o) having small 
L 3 j oc -norm of (uq, Vrfo)- Here ^ moc G^ 3 ) ' s ^ ne s P ace of uniformly locally L 3 -integrable functions. For any 
initial data (uo,do) with small ||(^o, V<io)ll £3(183): we show that there exists a unique, global solution to 
111. It which is smooth for t > and has monotone deceasing L 3 -energy for t > 0. 
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1. Introduction 



In this paper, we consider the Cauchy problem for the following hydrodynamic system modeling the flow 
of nematic liquid crystal materials in R 3 : for < T < oo and (u,P,d) : R 3 x [0, T) -> R 3 x K x S 2 , the 
system is given by 

' u t +u-Vu- vAu + VP = -XV ■ (Vd Vd), in R 3 x (0, T), 
V-u = 0, inR 3 x(0,T), 

(1 1) 

d t +u- Vd = 7 (Ad+ |Vd| 2 d), inR 3 x(0,T), 
(u, d) = (u , d ), on R 3 x {0}, 



for a given initial data (u ,d a ) : K 3 — > R 3 x S 2 with V • u = 0. Here it : R 3 — > R 3 represents the velocity 

field of the fluid, d : R 3 — > S 2 - the unit sphere in R 3 - is a unit vector field representing the macroscopic 

molecular orientation of the nematic liquid crystal material, P : R 3 — > R represents the pressure function. 

The constants z/, A, and 7 are positive constants that represent the viscosity of the fluid, the competition 

between kinetic and potential energy, and the microscopic elastic relaxation time for the molecular orientation 

field. V- denotes the divergence operator in R 3 , and Vd© Vd denotes the symmetric 3x3 matrix: 
OO 1 

O ; (Vd Vd) = (Vid, Vjd), 1 < *, j < 3. 



y 

Throughout this paper, we denote (v, w) or v ■ w as the inner product in R 3 for v, w £ R 3 . 

The system (jTTTJ) is a simplified version of the famous Ericksen-Leslie model for the hydrodynamics of 
nematic liquid crystals developed by Ericksen and Leslie during the period of 1958 through 1968 [6| fTol 14]. 
This system reduces to the Ossen-Frank model in the static theory of liquid crystals. It is a macroscopic 
continuum description of the time evolution of the materials under the influence of flow field u and the 
macroscopic description of the microscopic orientation field d of rod-like liquid crystals. The current form of 
system (jl.ll) was first proposed by Lin [17] back in the late 1980's. From the mathematical point of view, 
(jTTTJ) is a system coupling the non-homogeneous incompressible Navier-Stokes equation and the transported 
heat flow of harmonic maps to S 2 . Lin-Liu [19, 20 initiated the mathematical analysis of (JTTTJ) by considering 
its Ginzburg-Landau approximation or the so-called orientation with variable degrees in the terminology of 

Ericksen. Namely, the Dirichlet energy J -|Vd| 2 for d : R 3 — > S 2 is replaced by the Ginzburg-Landau 

f 1 1 

energy / -|Vd| 2 H r(l - \d\ 2 ) 2 (e > 0) for d : R 3 ->• R 3 . Hence (JTTTJb is replaced by 

J 2 \t l 

d t d + u-Vd = -f(Ad+ \(l-\d\ 2 )d). (1.2) 

Lin-Liu proved in [191 120] (i) the existence of a unique, global smooth solution in dimension two and in 
dimension three under large viscosity i>; and (ii) the existence of suitable weak solutions and their partial 
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regularity in dimension three, analogous to the celebrated regularity theorem by Caffarclli-Kohn-Nirenberg 
[3] for the three-dimensional incompressible Navier-Stokes equation. 

As already pointed out by [19, 20 , it is a very challenging problem to study the convergence of solutions 
(u e ,P e ,d e ) to (|l.ip i- (jl.ll) 9- (|1.2l) when e J, 0. In particular, the existence of global Leray-Hopf type weak 
solutions to the initial and boundary value problem of (11.11) has only been established recently by Lin-Lin- 
Wang [21] in dimension two, see also Hong [9] and Xu-Zhang [29] and Hong-Xin 12 for related works. 

Because of the super-critical nonlinear term V ■ (Vd V) in (jl.ljl i . it has been an outstanding open 
problem whether there exists a global Leray-Hopf type weak solution to (|l.ip in R 3 for any initial data 
(u ,d ) £ L 2 (R 3 ,R 3 ) x H Al ^ 2 (R 3 ,S' 2 ) with V • u = 0. It is standard that in R 3 the local existence of a 
unique, strong solution to can be obtained for any initial data u € VF S ' 2 (R 3 ) and d Q £ W S+1 ' 2 (M. 3 , S 2 ) 
for s > 3 with V ■ uq = 0, see for example 28 . A blow-up criterion for local strong solutions to (|l.lj) . similar 
to the Beale-Kato-Majda criterion for the Navier-Stokes equation (see [I]), was obtained by Huang- Wang 
[IT] . For small initial data in certain Besov spaces, Li- Wang [33] obtained the global existence of strong 
solutions to (jl.lj) . We would like to mention that Wang [37J has recently obtained the global (or local) 
wcll-posedness of (|1.1[) for initial data (ito, do) belonging to possibly the largest space BMCT 1 x BMO with 
V • Mo = 0, which is a invariant space under parabolic scaling associated with (II. ip . with small norms. 

In this paper, we are mainly interested in the local well-posedness of for any initial data (uo,do) 
such that (u , Vdo) € L 3 loc (R 3 ). Henceforth L 3 j qc (R 3 ) denotes the space of uniformly locally L 3 -integrable 
functions. It turns out that L^ qc (R 3 ) is also invariant under parabolic scaling associated with 

Now we give the definition of (1BL 3 ). The readers can consult the monograph by Lemarie-Rieusset 

[IB] for applications of the space Z/^ (R 3 ) to the Navier-Stokes equation. 

Definition 1.1. A function f E L 3 oc (R 3 ) belongs to the space £ 3 i oc (K 3 ) consisting of uniformly locally 
L 3 -integrable functions, if there exists < R < +oo such that 

ll/IU (K3) := sup ( f |/| 3 )' < +oo. (1.3) 

It is clear that 

. L 3 (R 3 ) C L 3 u1qc (R 3 ). 

• If / G L 3 j Qc (R 3 ), then \\f\\L 3 R (v^) is finite for any < R < +oo. For any two < i?i < R 2 < oo, it 
holds 

ll/ll^s) < \\f]\ L%2m < (g) , V / E L 3 ulQC (R 3 ). (1.4) 

' L uloc( R3 ) c D BMOfl 1 ^ 3 ) ( see US or [37]). Moreover, for any < R < oo, it holds 

0<R<oo 

[/]bmo-( K 3) < \\f\\ Ll(m , V / E L 3 ulQC (R 3 ). (1.5) 

Throughout this paper, we write A < B if there exists a universal constant C > such that A < CB. 
Here are a few more notations and conventions that we will use through this paper. For two matrices M, N 
of order 3, we use M : N — M^N 1 ^ to denote their scalar product. For two vectors u, v E R 3 , we 

l<ij<3 

let u ® v denote their tensor product: (u ® = u l v J , 1 < i,j ' < 3. For < s < +oo and 1 < p < oo, we 
denote by W S ' P (R 3 ) and W S ' P (R 3 as the Sobolev space and the homogeneous Sobolev spaces respectively. 
For < a < b < +oo, denote 

C fc °°(K 3 x [a, b}) = p| {/ £ C m (R 3 x [a, b}) : ||/||c»(r»x[o,6]) < +^}- 

m>0 

^([a^L^R 3 )) = {/ E L^([a,b],L 3 (R 3 ))}, 

and 

C°([a,6],L 3 loc (R 3 )) = {/ E C((a, b], L 3 (R 3 )) n L°°([a, b], i 3 (R 3 )) : as 1 1 0,/(t) -> f(a) in if oc (» 3 )}- 
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Repeated indices are summed unless specificized otherwise. Upper indices denote components and lower 
indices denote derivatives. 

Now we state our main theorem. 

Theorem 1.2. There exist e > and r > such that if u a : M 3 — > R 3 7 with V • uq = 0, and do : R 3 — > S 2 
satisfies (do — eo) € i 3 (M 3 ) /or some eo € 5* 2 , and 

|||(«o,Vdo)|||z»(ji») := sup ( / K| 3 + |Vd | 3 ] <e (1.6) 

for some < R < oo, then there exist To > to-/? 2 and a unique solution (u, d) : R 3 x [0, To) -> K 3 x M x S 2 

o/ sucft £/ia£ the following properties hold: 

(i) For 1 1 0, (u(t),d(t)) -> (u , d ) and Vd(i) -> Vd in £? oc (R 3 ). 

(u,d)e f) Cr(R 3 x[^„-i],M 3 x5 2 ), (u,Vd)e f| C°([0,T'],L 3 loc (K 3 )). 

0<(5<T 0<T'<T 

|||(u(t), Vd{t))\\\ L ce ([0tToR2]tL 3 Bm) < Ce . (1.7) 
(iwj //To < +oo is t/ie maximum time interval then it must hold 

UmBup|||(«(i),Vd(t))||| i s ( Hs ) >eo, V < r < oo. (1.8) 



The ideas to prove Theorem 11.21 are motivated by those employed by 21J. There are five main ingredients, 
which include 

• approximate (uo,do) by smooth (tip, dp) (see Lemma [5.11 below) and obtain < Tk < +oo and a 
sequence of smooth solutions (u k ,P k ,d k ) of (jl.ip in R 3 x [0, T k ], under the initial data (uq, d^); 

• utilizing the local L 3 -energy inequality (|3.ip . obtain uniform lower bounds of 

• apply the eo-regularity Theorem 14.41 to obtain a priori derivative estimates of (u k , d k ) and then take 
limit to obtain the local existence of L 3 loc -solutions to (|1.1|1 ; 

• apply Theorem 14. 41 again to characterize the finite maximal time interval; and 

• adapt the proof of [57] to show the uniqueness. 

For a solution (u,P,d) to (|1.1[) , denote its L 3 -energy by 

£7 3 («,Vd)(t)= I (K*)| 3 + |Vd(i)| 3 ), t>0. 

JR 3 

Concerning the global well-posedness of ()1.1|) , we have 

Theorem 1.3. There exists an e > such that if {u Q ,d ) G L 3 (R 3 ,R 3 ) x W 1 ' 3 (M? , S 2 ) , wzt/i V • u = 0, 
satisfies 

E 3 (u 0> Vd ) <e 3 , (1.9) 
i/ien i/iere exists a unique global solution (u, d) : R 3 x [0, oo) I 3 x R x S 2 o/ (H2]) such that (u, d) G 
C°°(R 3 x (0,+oo)) n C([0,oo),L 3 (R 3 ) x W^ 3 (R 3 )), E 3 (u,Wd)(t) is monotone decreasing for t > 0, and 

l|V"Xt)|| ioo(K3) + ||V m+1 d(t)|| LOO(R3) < V 1 > 0, m > 0. (1.10) 

We mention here that the first conclusion of Theorem 11.31 has been proven by [5] , which is based on 
refinement of the argument by Wang [27]. Since the exact values of v, A, 7 don't play a role in this paper, 
we henceforth assume 

v = A = 7 = 1. 

The paper is written as follows. In §2, we derive an inequality for the global L 3 -energy of smooth solutions 
of (jl.ip . In §3, we derive an inequality for the local L 3 -energy of smooth solutions of (jl.lj) and prove Theorem 
11.31 In §4, we will prove an eo-regularity for suitable weak solutions to {JTTTj). In particular, a priori derivative 
estimates hold for smooth solutions to (|1.1[) under a smallness condition. In §5, we will prove Theorem II .21 
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2. Inequality on the global L 3 -energy and proof of Theorem 11.31 
In this section, we will derive an inequality for the L 3 -energy £3 (it, Vd) (i) for any smooth solution 



(u,d) 



[0,T] 



x § 2 , for < T < 00, of the system ([1.11) for nematic liquid crystals. 



Lemma 2.1. There exists C > such that for < T < 00 if (u,d) e C°°(R 3 x [0,T),M 3 x S 2 ) n 
C([0, T), L 3 (M 3 ) x W /1 ' 3 (R 3 )) anrfPeL°°([0,T),Li(K 3 )) solves fT7]) . i/ien « Mds 



(2.1) 



- /_(| u | 3 + |vd| 3 )+ 1 - c|| u ||i3 (B3) /_H|Vi»| 



1 - C(||U|| L 3 (M 3) + |M|l3( R 3)|!W|| L 3 (K 3) + ||Vd||i3 (R 3)) / \V d\\V ' d\' < 



Proof. Taking spatial derivatives of (jl.ip ^, multiplying the resulting equation by |Vd|Vd, and integrating 
over R 3 , we have 

— f I|Vd| 3 =/ V(Ad) : |Vd|Vd - [ V(u ■ Vd) : \Vd\Vd - [ V(|Vd| 2 d) : |Vd|Vd. (2.2) 

dt Jm3 3 J R 3 J R 3 J R 3 

For terms on the right hand side of ([2.21 by integration by parts we have 



V(Ad) : |Vd|Vd 



! n{|Vd|>o} 



V 2 d : V(|Vd|Vd) 



< 



(|Vd||V 2 d| 2 
|Vd||V"d| 2 



>n{\Vd\>0} 

72 7|2 



|V 2 d- Vdp 
|Vd| 



/ V(w • Vd) : |Vd|Vd = - / (u-Vd)-((V|Vd|)-Vd+|Vd|Ad) 

Jr 3 Jr3 

< / |u||Vd| 2 |V 2 d|, 

and, using \d\ = 1, 

/ V(|Vd| 2 d) : |Vd|Vd = / (V|Vd| 2 ) • |Vd|V(^-) + / |Vd| 2 Vd : |Vd| 2 Vd 

Jr 3 Jr 3 2 j R 3 

= f |Vd| 5 . 

JR 3 

Putting these estimates into ()2.2[) yields 



±[ |Vd| 3 +/' |V(|Vd|l)| 2 < / |Vd| 5 + M|Vd| 2 |V 2 d|, 

Jr 3 Jr 3 Jr 3 



where have used the following variant of the Kato inequality 

|V|Vd|* 



^|Vd|^|V|Vd||<^|Vd|5|V 2 d|. 



Observe that by the Sobolev inequality and the Kato inequality above, we have 



|Vd| 9 = / (|Vd| 5 ) 6 < / |V|Vd|i 
Jr 3 \Jr 3 

Hence, by the Holder inequality and ([2.4)1 . we have 



3 

12 1 < 



|Vd||V 2 d| : 



|Vd|| 5 L5(R3) < ||Vd||i3 (R 3 ) ||Vd|| 3 9(R 3 ) < (J |Vd| 3 ) ' (J |Vd||V 2 d| 



(2.3) 



(2.4) 
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For the second term on the right-hand side of (|2.3|) . by the Holder inequality and (|2.4[) we have 

M|Vd| 2 |V 2 d| < \\u\\ L3(m \\\Vd\i\\ L e m \\\Vd\^\V 2 d\\\ L 2 (R3) 

< ||u||L3 (R 3 ) ||V|Vd|^|| L 2 (R 3 ) |||Vd|^|V 2 d||| L 2 (R 3 ) 

< \\u\\ LSm \\\Vd\i\V 2 d\\\ 2 L2m . 
Inserting these two estimates into (12.3[) yields 

JtS 3 |W|3 + I 1 " C (" W II' 3 (K 3 ) + NU'(*»))] ^ 3 |Vd||V 2 d| 2 < 0. (2.5) 
Next we estimate the L 3 -norm of u. Multiplying by \u\u and integrating over R 3 gives 



d r i, l3 



dt /in 3 3 



Au ■ \u\u - / (u ■ Vu) • \u\u - VP • \u\u - (V ■ (Vd Vd)) • |u|tt. 

7R3 JR3 JK3 

For the terms on the right hand side of (|2.6I) . by integration by parts we have 

(Au) ■!«!«=- / |Vu| 2 | u | + | u ||V|«|| 2 , 



f (u-Vu)-\u\u= [ u-v( l ^pj =0, 

JW 3 J R 3 \ / 

/ VP-\u\u = - Pu-V\u\+P\u\(V -u) = - Pu-V\u\, 

JR3 Jr3 Jr3 



and 

r (V ■ (VdQVd)) -\u\u= I (VdOVd) : V(|«|u) 



< 



(Vd© Vd) : V|u| ®u + |u|(Vd© Vd) : Vu 
|Vd| 2 |u||Vu|. 



Substituting these estimates into (|2.6j) . we obtain 

/ ' 3 4- / l7/.IIV7/J 2 |7/.l < 



(2.6) 



u\\Vu\ 2 \u\< / |P||u||V|u||+ / |Vd|>||Vu|. (2.7) 



Using the Kato inequality |V|«|| < |V«|, the Cauchy inequality and the Holder inequality in (|2.7I) . we obtain 
±[ |u| 3 +/ M|V U | 2 <c/ \ u \(\P\ 2 + \Vdf) + U \u\\Vu\ 2 

at JR 3 JR3 JR3 i J R 3 

< c(II^IIl3( H 3) + l|v<f||i. (R 3 ) )|MU3 CR 3 ) + \j \u\\vu\ 2 . 

Therefore we get 

±[ \u\ 3 +f H|Vu| 2 <(||P||i, (RS) + ||Vd||i. (R . ) )||u|| i » (R .). (2.8) 

"t JR3 JR3 

We need to estimate ||P||^3(R3). To do so, we take divergence of (|l.llh to obtain 

-AP = V • V • (u®u + VdQ Vd). (2.9) 

Set 

g jk := v?u k + Vjd ■ V fe d, 1 < j, k < 3. 

Then we have 

P = A-^V^'*) - -RjRfcO^*). (2.10) 
Henceforth Rj = (— A)~ 5 V.,- denotes the j t,l -Riesz transform on M 3 for 1 < j < 3. 
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Since It, : L q (R 3 ) -> L q (R 3 ) is bounded for 1 < q < oo (see Stein [25]), we have 

||P|U3 (R3) = IIR.RfcCg^)^.,^) < ||^ fe |U3 (K3) < ||u||| 6(K3) + ||Vrf|||a (K3) . (2.11) 
Inserting (pmj) into (|2J| yields 

4/ M 3 +/" |«||Vu| 2 <(||«||i. (R , ) + ||Vd||i« (R 3 ) )||«||L.(R3). (2.12) 
«* Jr 3 Jr 3 

Using the Holder inequality, the Sobolev inequality, and |V|it|'| < |Vu||u|', we have 

\\u\\h'(R 3 ) < l|w||L3(R 3 )ll u llio(E3) < ||w|U3(R3)||V|n|5||| 2(a3) < \\u\\ L s iM3) J ^ |u||V|u|| 2 . 

Similarly we have 

l|Vd||ia (R 3) < ||Vd|| L 3 (R3) f |Vd||V 2 d| 2 . 

Jr 3 

Substituting these two estimates into (I2.12[) , we obtain 

±[ \u\ 3 +[ \u\\Vu\ 2 < \\uf L3m [ |«||Vt*| 2 + ||«||isc»»)ll Vd lU-(» a ) / |Vd||V 2 d| 2 . (2.13) 
at Jr 3 Jr 3 Jr 3 Jr 3 

Combining and ([235]) yields pTTjl . □ 

Corollary 2.2. T/iere exists e > suc/i that for < T < oo ; if (it,d) £ C°°(JR 3 x [0,T),R 3 x S 2 ) D 
L°°([0,T),L 3 (R 3 ) x W 1,3 (R 3 )) is a solution to {OP satisfying 

E 3 {u ,Vd )<e 3 , (2.14) 

t/ien E 3 (u(t), V<i(t)) is monotone decreasing for < t < T . 

Proof. Denote 

E 3 {t) := E 3 (u(t),Vd(t)), t> 0. 

Let T max e [0, T) be defined by 

T max = max {t e [0, T) : E 3 {s) < 2e 3 , V0<s<<}. 

By continuity and (|2.14l) . we have that < T max < T, and 

E 3 (t) < 2e 3 , 0<t< T max , E 3 (T max ) = 2e 3 . (2.15) 
Suppose T max < T. Choose eo > so small that 

1 - Cel > and 1 - C*(e + 2eg) > 0. 

Then (f2TT5|) and pTTjl imply that 

< ^ 8 (*) + [1 - C6 2 ] / u\\Vu\ 2 + [1 - C*(e + 2eg)] f \Vd\\V 2 d\ 2 < 

at at j R 3 Jr3 

holds for < t < T max . Hence E 3 (t) is decreasing in [0,T max ] and 

E 3 {T m ^) < E 3 (0) <e 3 < 2e 3 . 
This contradicts the definition of T max . Thus T max = T and E 3 (t) is monotone decreasing in [0,T). □ 

Proof of Theorem [Ell Since C°°(M 3 , S 2 ) is dense in W^ 3 (M 3 , S* 2 ) (see [24]), it is not hard to show that 
there exist { (itg , d§ ) } C C°° (R 3 , R 3 ) x C°° (M 3 , S 2 ) such that 

V ■ u k = in R 3 , lim (||ttg - Mollis™) + ||V(d§ - do)||LS(RS)) = 0. 

Consider the system (jl.ljl under the initial condition (u,<i)| t= o = (ug,dg)- It is standard that there exist 
T k > and smooth solutions (u k ,d k ) £ C*°°(R 3 x [0,T fc ],M 3 x S 2 ) H C([0, T k ], L 3 (R 3 ) x H /1 ' 3 (R 3 )) to (10) . 

Since E 3 (uq, Vdo) < e 3 ,, we may assume that S3 (tig, Vd§) < 2e§ for all fc > 1. Hence by Corollary |2.21 we 
conclude that 

sup E 3 (u k {t),S7d k (t)) < E 3 {ul,Vdl) < 2e 3 , V k > 1. 

0<*<T fc 
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For the corresponding pressure functions P k , since 

Ap k = _v • V • (u fe <g) u k + Vd k Vd k ) in R 3 , 

we have 

sup ||P fe || 3 < sup (h fc ||i3 (R3) + ||Vd fe ||i3 (R 3 ) ) <Ce 2 . 

0<t<T k L (K > 0<t<T k 

Let T k be the maximal time interval for (u k ,d k ). If < < +00, then by Theorem 14.41 in §4 below we 
conclude that (u k ,d k ) G C b °°(R 3 x [0,T fe ],R 3 x S 2 ). Hence (u k {T k ),d k (T k )) G C°°(R 3 ,R 3 x S 2 )nL 3 (K 3 ) x 
W h3 (R 3 ), and 

E 3 ((u k (T k ),Vd k {T k ))) < 2e 3 

so that we can extend the smooth solutions (u k , d k ) beyond the time T k , which would contradict the maximal- 
ly of T k . Therefore T k = 00 and the smooth solution (u k ,d k ) exists globally. Moreover, E%(u k {t), Vd k (t)) 
is monotone decreasing and less than 2eg. By Theorem 14.41 we have the derivative estimates: 

HV ro u fc (*)|| ioo(R3) + ||V™ +1 4(*)|U (R3) < V t > 0, m > 1. (2.16) 

After taking possible subsequences, we may assume that there exists (u,d) G C°°(R 3 x (0, +00), R 3 x S 2 ) n 
C([0, +oo),i 3 (R 3 ) n VK 1 ' 3 (R 3 )) such that as k -> 00, 

(1) (u k ,d k ) -> (u,d) in C£ C (R 3 x (0,+oo)) for any m > 1. 

(2) (u fe , Vd fe ) ->• (u,Vd) weak* in L°°([0, +00), L 3 (R 3 )). 

Thus (u,d) G C"°°(R 3 x (0,+oo),R 3 x S 2 ) solves fTTjh. (fTT]) ? , and (HUDa, and the estimate PTTU]) holds. 
Using the equation (jl.ip . we can get that for any < T < +00, 



sup 

fe>i 



{d t u k , d t d k ) 3 < C(T) < +00. 

L2([0,r),vr _1 '2 (R 3 )) 



This implies that (u, d) G C([0, T], i 3 (R 3 ) x Ty 1 ' 3 (R 3 )) and (u,d)| t=0 = (u ,d ). Applying Corollary 
again, we conclude that E%(u(t), VcZ(t)) is monotone decreasing for £ > 0. The part of uniqueness can be 
proved as in the step 6 of the proof of Theorem ll.2l in §5, which is omitted here. The proof is complete. □ 

We would like to mention applications of Theorem II .31 to the heat flow of harmonic maps and the Navier- 
Stokes equation. 

1) If u = 0, then reduces to the heat flow of harmonic maps to S 2 for d : R 3 x (0, +00) — s- S 2 : 

\d t d = Ad+ \Wd\ 2 d in R 3 x (0, +00) 

\ d = d Q on R 3 x {0}. ( ' ' 

2) If d is a constant unit vector, then (|l.l[h and (| 1 . If) -? reduce to the Navier-Stokes equation: 

'dtU + u-Vu- Au + VP = in R 3 x (0, +00) 

V-m = inR 3 x(0,+oo) (2.18) 
u = u Q on R 3 x {0}. 



The following properties follow directly from Theorem 11.31 We would like to point out the observation of 
monotone decreasing property of the L 3 -energy seems new. 

Remark 2.3. 1) There exists eo > such that if do : R 3 — > S 2 satisfies / |Vc?o| 3 < then there is a 
unique global solution d : R 3 x [0, +00) -> S 2 of (j2~T7)) such that d G C([0, +00), VK 1 ' 3 (R 3 , S 2 )) H C7°°(R 3 x 
(0, +00), S 2 ), and / |Vc?(t)| 3 is monotone decreasing for t > 0. 

2) There exists eo > such that if uq : R 3 — > R 3 , with V-i*o = 0, satisfies / |uo| 3 < £q, then there is a unique, 
global solution u : R 3 x [0, +00) -> R 3 of ([2~T5|) such that u G C([0, +00), L 3 (R 3 )) n C°°(R 3 x (0, +00), R 3 ), 

and / |u(i)| 3 is monotone decreasing for t > 0. 

Jr 3 
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3. Inequality of the local £ 3 -energy 

In this section, we will derive an inequality of the local L 3 -energy for smooth solutions (u, d) : M 3 x [0, T] —> 
M 3 x S 2 for < T < oo , of the system (jl.ip . More precisely, we have 

Lemma 3.1. There exists C > such that for < T < oo, if (u,d) € C°°(R 3 x [0,T),R 3 x S 2 ) H 
C([0,T),L 2 (R 3 ) x W h2 (R 3 )) is a smooth solution of the system iTO]) , then 

± [ (H 3 + |vd| 3 )0 2 + / (\v(\u\ 3 ^)\ 2 + \v(\vd\ a H)\ 2 ) 

at JR 3 J R 3 V / 

<C [ (M 3 + |Vd| 3 )|V0| 2 + CiT 2 sup ( f |u| 3 + |Vd| 3 ) (3.1) 

JR 3 y£R 3 \JB R (y) J 

2 

+ c([ \u\ 3 + \Vd\ 3 ) 3 f (\V(\u\h)\ 2 + |V(|Vd|i0)| 2 ) , 

\Jspt<j> / JR 3 v ' 

holds for any4>^ C(f(R 3 ), withO < <j> < 1, spt <f> = B R (x f\ for some R > and x Q € R 3 , and \V<j>\ < ABr 1 . 
Proof. We divide the proof into three steps. 

Step 1. Estimation of the local L 3 -energy of Vd Differentiating (jl.llh with respect to x, integrating against 
(f) 2 \Vd\Vd over R 3 , and applying integration by parts, we have 

d 



dt 



\Vd\ 3 cj) 2 + 3 f V 2 d: V(ct> 2 \Vd\Vd) 
Jr 3 



(3.2) 



<3 f \\7d\ 5 cj) 2 + 3 [ (u-Vd)-V-(0 2 |Vd|Vd), 

JR 3 JR 3 

where we have used \d\ = 1 and the following identity to obtain the first term on the right hand side: 

V(|Vd| 2 d) • \Vd\(Vd) = ^V(|W| 2 )|W|V(|d| 2 ) + \Vd\ 3 Vd ■ Vd = |W| 5 . 

For the second term on the left hand side of (|3.2j) . direct calculations using |V|V<i|| < |V 2 d| and the Holder 
inequality imply 

f V 2 d:\7(<j) 2 \Vd\Vd)= [ \Vd\\V 2 d\ 2 (j) 2 + [ (|Vd| 2 V|Vd| • V0 2 + | Vd| | V| Vd| | 2 ^> 2 ) 

JR 3 JR 3 JR3n{|V(i|>0} 

>~ / |Vd||V 2 d| 2 2 -C f |Vd| 3 |V0| 2 . 

A Jr 3 Jr 3 

For the second term on the right hand side of (|3.2j) . by the Cauchy inequality we have 
(u ■ Vd) ■ V • (^ 2 |Vd|Vd) < 2 f \u\\\7d\ 2 \\7 2 d\<f> 2 + M|Vd| 3 </>|Vc£| 



<\[ \vd\\v 2 dU 2 + c(f nAUf |Vd|V 6 

8 Jr 3 \Jspt<j> J \Jr 3 

i i 
\ 3 / r \ e / r 

3 II |V7J|9 /6 \ I / |V7J|3|V7 j,|2 



C ( / \u\ 3 ) ( / |Vd| y ^ ) ( / |Vd| 3 |V0| 

/ Spt0 



< - 

~ 8 



I f |Vd||V 2 d|V + C f |Vd| 3 |V^ 

° Jr 3 Jr 3 



C{ I \u\ 3 ) ( / |Vd| 9 - 6 

' Spt0 



By the Holder inequality and the Sobolev inequality, we have 



|Vd|V 6 ) </ |V(|Vd|M| 2 , / |Vd| 5 2 < ( / |Vd| 3 ) / |V(|Vd|^)| 2 

' spt0 



^Here spt</> denotes the support of * 
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Putting these estimates into (|3.2p yields 

d_ 

dt 



r|Vd| J + / |V 2 d| 2 |Vd|0 2 



< 



|W| 3 |V0| 2 



|Vd| ; 



Spt0 



V(|Vd|5</>) 



(3.3) 



Step 2. Estimation of the local L 3 -energy of u. Multiplying (fOJi by 2 |u|u and integrating over R 3 yields 



dr 

dt J R 3 



< 



|m| 3 ^ + 3 / Vu-V(^\u\u) 



\Vd\\V z d\\u\ z (t, z + / |Vu||u|V + / |P-c||V(^|«|u)| 



where c <E K is a constant to be chosen later. 

By the Cauchy inequality, the Holder inequality, and the Sobolev inequality, we have 

\7u-\7(tf\u\u) > 1 / M|V U |V~4 / M 3 |V0| 2 , 



(3.4) 



|Vu||u| 3 2 < - / |w||Vu|V + C / |u| 



Spt0 



V(|u|*0) 



and 



/ |Vd||V 2 d|| ? 



"|V<§/ |Vd||V 2 d| 2 2 + C / |Vd||u|V 



< - / |Vd||V 2 rf| 2 2 + C / |Vd| 



Spt0 



spt0 



|u|V 



For the last term on the right hand side of (|3 .4[) we have 

\p - C \\v ■ (\u\utf)\ <l f u\\vu^ + c I |.r-H' J l4> J • / M- ! |v,.| 



Putting these inequalities into (|3.4p we obtain 



5' |u 



3^2 



|w||Vu|V 



1 



<C^J M | 3 |V0| 2 + ^ / |V ( /||V- ( /|-m- -C I \P ~ c\ 2 \u\<f> 2 



Combining 



spt</> 

with (|3.5p yields 
d 



|Vd| ; 



(|v(H^)| 2 + |v(|vd|^)| 2 ). 



(|u| 3 + |Vd| 3 )0 2 + / (|V(M^)| 2 + |V(|Vd|l0)|- 

./1B3 \ 



< c / (M 3 + |vd| 3 )|V0| 2 + c / M|p - c | 2 ^ 



+ c 



spt<f) 



u\ 3 + \Vd\A 3 f (|V(|u|*0)| a + |V(|Vd|^)| : 



Step 3. Estimation of the pressure function P. By the Holder inequality, we have 



H|P- C |V< \u\ 



Spt0 



|P-c 



3 j.3 



(3.5) 



(3.6) 



We see that ()3.1|) follows from (|3.6p and the estimate (|3.7p of Lemma I3~^l below. The proof is complete. □ 
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Lemma 3.2. Under the same assumptions as in Lemma \3.1\ assume that (j> € Co°(R 3 ) satisfies < cj) < 1, 
spt = Br(xq) for some xq £ R 3 , and \V(f>\ < 2R . Then there exists C > such that for any t e (0,T) 
there is c(t) el so that the following estimate holds 



\P(t)-c(t)\ 3 ^) <C( \u(t)\ 3 + \Vd(t)\ 



spt<p 

CR- 1 sup 



(\V(\u(t)\H)\ 2 + \V(\Vd(t)\H)f 



\u(t)\ 3 + \Vd(t)f 



V m 3 \JB R (y) 

Proof. For simplicity, we write (w, P, d) and c for (u(t), P(t),d(t)) and c(t) respectively. Since 



AP = V%(gJ k ), g jk := u>u k + V 3 d ■ V,,,/. 



we have 



P = -R i R fc (^' fc ) 
where Rj is the j-th Riesz transform on R 3 . Hence we have 

(P-c)0 = -R,R fc (^' fc )0-c0 

= -R j R k (g ik <l>) - [4>,IL j K k ]{ff lk ) ~ c(f> 
where [(p, RjRfc] is the commutator between <p and R^Rfc given by 

[0, 11,-RfcK/) = • R,R fc (/) - R,R fe (/0), / e C °°(R 3 ). 
We now estimate [0, R^-R^] (g J ) as follows. 

[^RjRfc] ( 5 jfe )(z) 

= 0(z)RjR fc (<? ifc )(z) - RjRfc^'V)^) 



0(x 



J'Vt 1 -,,^ ., r ( x i - y j )(x k ~ y k ) 



J»3 \x-y\ 5 

(</>(x) - 0(y))(a;J - y 3 ){x k - y*) - fc 

I is y vi//* 



For any x £ spt = Br(xq), we have 

[^RjRfe] (^' fc )(x) + c#c) 



B2i?Oo) 



[5 9 {y)dy + c<p(x) 



\x - y\ 
{ X i - yi)(x k - y k ) nik 



I(x)+II(x). 



\x - y\" 



g J *{y)dy + c 



For I(x), we have that 



\4>{x) - (j>{y)\\x3 - yi\\x k - y k \ jk 
B 2R (x ) \x-y\ 5 

XB 2RM (y)(\u\ 2 + \Vdf)(y) 



(3.7) 



(3.8) 



< CR , 

\x - y\ z 

= CR-% ((\u\ 2 + \Vd\ 2 ) XB2R(xo) ) (x), 

where Xb 2R (x ) is the characteristic function of B>2r(xo), and Ii is the Riesz potential on R 3 of order 1 given 
by 



Uf)(x) 



\f(y)\ 



, x £ R , V / G L (R ) 



\x - y\ 2 '' 

Recall that by the Hardy-Littlewood-Sobolev inequality, Ij : Li (R 3 ) — > L 3 (R 3 ) satisfies 

I|Ii(/)IIl3( R 3) < ll/ll,j (R3) - 



(3.9) 
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Hence we have 



,) < R- 1 Hi! ((|u| 2 + \Vd\ 2 ) X B 2R (x )) 

< R-'WQul 2 + \Vd\ 2 ) XB2R{xo) \\ L s 



f \u\ 3 + \Vdf 

J B 2R (xo) 



(3.10) 



< R- 1 sup [ f 

i/6l 3 \ J E 



+ |Vd| 2 



To estimate 77, choose 

Note that 

Therefore we have 

\II(x)\ = 

< 



yeS. 3 \JB„(y) 

(x Q -y) j (x a -y) k jk 



i(x ) \xo - y\ 5 



g Jk (y)- 



3,k 



< +0O. 



<K*) [ 


3 \b 2R (x ) ( 







R\<fi(x)\ 



Ir 3 \b 2R ( X0 ) \x~v\ 4 

where we have used the following inequality (see |25j): 

(x j -y j ){x k -y k ) (x - y) j (x - y) k < \x - x\ 



{xi ~yi)(x k -y k ) _ (x - y)i (x Q - y) k \ ^ fe 
\x-y\ 5 \x -y\ 5 

' (\u\ 2 + \Vd\ 2 )(y), 



\x - y\* 



Thus we have 



l^o - y\ 5 

\II\{x)<r( 
Jr 



\xo-y\ 4 
1 



, for x e B R (x ) and y e M 3 \ B 2R (x ). 



< 



w\b 3R (xo) \ x ° - y\ 4 



(H 2 + |Vd| 2 )(y) 

(M 2 + |W| 2 ) 



B (k+1)R (x )\B kR (x ) 



< 



R 3 



< R~ 2 



00 

sup / (M 2 + |V<2| 2 

,fc=2 J veR 3 JB R (y) 



sup 

v m 3 \ JB R (y) 



\u\ 3 + |Vd| E 



Integrating 7/ over Br(x ) we get 

||77|U3 (R 3) < Ti" 1 sup ( / M 3 + |W| 3 

y£R 3 \JBr(v) 



Additionally, we have 



|R J R fe (^ fe 0)|| i3(K 3 ) 



< 



< 



< 



y£R 3 \JB R (y) 

(| u | 6 + |Vd|V 
\u\ 3 + |W| 3 



(3.11) 



Spt0 



+ |Vd| : 



sptcp 



(|u| 9 + |Vd|V 6 
(|V(|u|l0)| 2 + |V(|Vd|^)P 



(3.12) 



Combining the estimates (|3.10[) and p. lip with (|3.12l) yields (|3.7[) . This completes the proof of Lemma 
O □ 
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4. Regularity of Suitable Weak Solutions 

In this section, we will derive a priori estimates for smooth solutions to the system Ijl.ljl under a smallness 
condition for the L 3 -norm of (it,Vd). Since the method is flexible enough, it also yields the smoothness 
for a subclass of suitable weak solutions to (jl.ip . We present the result in the context of suitable weak 
solutions to (|l.ip . The notion of suitable weak solutions was first introduced by Caffarclli-Kohn-Nirenberg 
[3] in the context of incompressible Navier-Stokes equations. Here we adapt this notion to ([1.1)1 . similar to 
the definition given by Lin )18) on the Navier-Stokes equation. 

Let < T < oo and f2 C K 3 be a bounded smooth domain. 

Definition 4.1. A triple of functions (u, P, d) : fl x (0, T) — > M 3 x K x S 2 is called a suitable weak solution 
to the system in x (0, T) if the following properties hold: 

(1) u £ LfL 2 , n L 2 Hl(Vt x (0,T)), PeLi(fix (0,T)) and d e L 2 H 2 x {n x (0,T)); 

(2) (u, P, d) satisfies the system in the sense of distributions; and 

(3) (u,P,d) satisfies the local energy inequality (|4.1[) . 

Now we would like to point out that the class of smooth solutions belongs to the class of suitable weak 
solutions to the system (11.11) . Let I 3 denote the identity matrix of order 3. 

Lemma 4.2. Suppose that (u,d) £ C°°(n x (0,T),M 3 x M x S 2 ) is a solution of IP[) infix (0,T). Then 
for any nonnegative cf) S C^°(fl x (0,T)), it holds that 



2/ (\Vu\ 2 + \Ad+\Vd\ 2 d\ 2 )(b< f (M 2 + |W| 2 )(0 t + A0) 

- [ {\u\ 2 + |Vd| 2 + 2P)u ■ V<j> 
2 f (Vd0 Vd- |Vd| 2 I 3 ) : V 2 



(4.1) 



<(0,T) 

2 1 Vd Vd : u V</>. 

/qx(o,t) 



Proof. Multiplying (|l.l[) i by it0 and integrating the resulting equation over fix (0,T) yields 



UfU(j)+ (u ■ Vu) ■ U(f>— Au ■ u(f> + VP ■ itij 

f2x(0,T) Jflx(0,T) Jflx(0,T) Jnx(0,T) 



/ Vd0 Vd : V(u<£). 



(4.2) 



<(0/ 

Applying integration by parts, the terms on the left hand side of ([4.2)1 can be estimated by 



/ u t ■ 

Jnx(o,T) 


' U<f. 


—f 


<(0,T) 




1 


/ (u ■ Vu) 

/ox(0,T) 


■ u4 


—J 


<(0,T) 


X | |2 




/ Au 

Jnx(o,T) 


■ u4 




<(0,T) 


|Vu| 2 0- 


Jnx(o,T) 


I VP 

jQx(0,T) 


■ u4 


--I 


<(0,T) 


P(u • V. 


i>)- 



For the term on the right hand side of (I4.2[) . we have 

/ Vd0 Vd : V • (u<f>) = / Vd Vd : [(Vu)<£ + u V<; 

Jnx(O.T) JQx(O.T) 
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Putting these identities into (|4.2p yields 

f - l -\u\ 2 ^ t + M)-{ l -\u\ 2 + P){u-Vcj>)+ f |Vu| 2 </> 

Jnx(o.T) z z Jnx(o,r) 

(Vd Vd) : [^Vu + u <g> V0] . 



'fix(0,T) 

Differentiating with respect to x and integrating against (Vd)<^>, we have 



/flx(O.T) Jnx(O.T) Jfix(0,T) 

For the first term on the left hand side of (|4.4j) . we have 



f (Vd) t : (Vd)0 = - f ~|Vd| 



2 0f 



'f!x(0,T) Jnx(O.T) 

Using (JTTTJ) 2 , we can simplify the second term on the left hand side of (|4.4[) into 

V(it • Vd) : (Vd)i/) = / "adj • d a (f) + / v? dj a ■ d a (f> 

Ox(0,T) Jnx(O.T) J!ix(0,T) 



/ Vu:Vd0Vd0+/ • V(|Vd| 5 

Jnx(o.T) Jnxio.T) 2 



'nx(o.T) ^nx(o.T) 

Vu:Vd0Vd0-/" -(u • V0)|Vd| 2 . 

ox(o,T) Jnx(ox) ^ 



For the term on the right hand side of (|4.4j) . differentiating |d| = 1 gives 

Vd-d = and Ad ■ d + | Vd| 2 = 0. 
Thus, by integration by parts we have 

/ V [Ad+ |Vd| 2 d] • Vd</> = - / [Ad+ |Vd| 2 d] • [Ad<j) + Vd • V<t>] 

Jqx(o.t) Jnxio.T) 



[ |Ad+ |Vd| 2 d| 2 </>- / Ad-(Vd-V^). 

Jnxto.T) Jnxio.T) 



(4.3) 



f (Vd) t : (Vd)<£ + f V(u • Vd) : (Vd)0 = f V [Ad + | Vd| 2 d] : (Vd)</> (4.4) 

Jnx(o.T) Jnxio.T) Jnx(o.T) 



'fix(O.T) JOx(0,T) 

f 

'fix(0,T) JOx(0,T) 

By integration by parts we have 

- ( Ad(Vd-V0)=/" (Vd Vd) : V 2 </> - / i|Vd| 2 A0 
Jnx(o,T) Jnx(o.T) Jnx(o.T) 2 

= / (Vd0 Vd- |Vd| 2 I 3 ) : V 2 </>+ / -|Vd| 2 A0. 

Jnx(o,T) Jnx(o,T) 2 

Inserting these identities into (|4.4[) yields 

/ [-^|Vd| 2 (0 t + Acj)) - i|Vd| 2 (w- + I Vw:Vd0Vd(/> 

Jnx(o,T) 2 2 Jnx(o.T) 

= [ (Vd0 Vd- |Vd| 2 I 3 ) : V 2 - / |Ad + |Vd| 2 d| 2 0. 

Jnx(o.T) Jnx(a,T) 

Combining P~3")) with (|4"3|) yields (|3~T|) , □ 

Corollary 4.3. Suppose that (u, P, d) : i7 x (0, T) — > R 3 x M x S* 2 is a suitable weak solution of the system 
il.l]) in Q x (0,T). Then for any nonnegative <fr € C°°(f2 x (0,T)) and < t < T, it holds 



(4.5) 



f (M 2 + |Vd| 2 )0 + 2 f (|Vw| 2 + |Ad+|Vd| 2 d| 2 )0 
Jnx{t} Jnx(o,t) 

<[ (| M | 2 + |Vd| 2 )(0 t + A0)+ f (|u| 2 + |Vd| 2 + 2P)u- V</> (4.6) 

JClxtO.t) Jtix(O.t) 



'fix(O.t) Jf2x(0,i 

(Vd Vd - |Vd| 2 I 3 ) : V 2 </> + 2 / Vd Vd : u V<f>. 
nx(o.t) Jnx(o.t) 
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Proof. For e > 0, let r) e G C°°([0, t]) be such that < r] < 1, 77 = 1 in [0,t - 2e], and r? = in [t - e, *]. ([Ojl 
follows by first applying ()4. lj) . with </> replaced by r] e (t)4>(x,t), and then taking e — > 0. 

Let C(3) > denote the best Sobolev constant of M 3 : 

C(3) :=inf{ ll ^ llL2(R3) : ^ / € C °°(]R 3 )}, 

and 23(3) > denote the constant in the following W 2, 2 -estimate: 



□ 

(4.7) 
(4.8) 



||V 2 /IIl 2 ( Bi) < 23(3)|| A/|| L2(Bl) + C\\Vf\\ wi , a(BBi) , V / G W 2 ' 2 ^). 
For z = (x , t ) £t 3 x (0, T) and r > 0, denote 

B ro (x a ) = {x G K 3 : |x - ac j < »"o}, P-o^o) = B ra (x ) x (t - rg,t ]. 
Now we are ready to prove the following eo _re gularity theorem. 

Theorem 4.4. For any 5 > 0, t/iere exists eo > smc/i that (u, P,d) : il x (0, T) — > M 3 x R x S 2 is a suitable 
weak solution to and satisfies, for zq — (xo,to) € f2 x (0,T) and P ro (zo) C O X (0,T), 



-Pr (20) 



+ r„ 



Vd 



^ (z ) 



|Vd| 3 



<eo, 



< 



1-5 



L?LUP ro (z )) C(3)X>(3)' 



t/ien (u,d) G C°°(Pio (zo),R 3 x S* 2 ), and t/ie following estimate holds: 

||(tt,d)|| C m ( P (jSD )) < C(m,r ,e ), V to > 0. 



(4.9) 
(4.10) 
(4.11) 



The crucial ingredient to prove Theorem 14.41 is the following decay lemma, which is analogous to that of 
the Navier- Stokes equations by [18] and [7]. 

Lemma 4.5. For any S > 0, there exist e > and 8 G (0, 5) swc/i that if (u, P, d) : Q x (0, T)-^K 3 xlx S* 2 
is a suitable weak solution of 1 1.1]) . and satisfies, for zq — (xo,to) G SI x (0,T) and P ro (zo) C SI x (0, T), 
feot/i J^.ff[ ) and |^.t0| ), t/ien it /icdds t/iat 



(Vo)" 



< 



fe r (zo) 
2 

pi*V 



(0 O r O )' 2 
|Vd| 3 



Pe r (zo) 



|Vd| = 



(4.12) 



Proof. By the invariance of under translations and parabolic dilations, it suffices to consider the case 
that zq = (0,0) and ro = 1. We will prove the Lemma by contradiction. Suppose that the conclusion were 
false. Then there would exist So > such that for any 6 G (0, 1) there are a sequence of suitable weak 
solutions (ui,Pi,di) of (|1.1[) in Pi, that satisfy 



Pi 



Pi 



|Vd,| 3 
Pi 

1-do 



= -> 0, 



and 



|Vdi|| £f)£ , (ft) < c(3)2?(3) 



IPI« 



|Vd, 



1 

>2 



Pi 



1^1* 



Pi 



|Vdi| 



Pi 



(4.13) 
(4.14) 

(4.15) 
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Now we define the blow-up sequence («t, Qj, e t ) : A -> R 3 x 1 x I 3 by 

Vi(z) = , Qi{z) = , e l {z) = , z e -Pi, 



€i €i 



l r 

where (e?;)i = j—r / di is the average of di over P\. 
I"il J Pi 

Then (vi, Qi, et) satisfy the following equations in Pi: 

' d t v, - A Vl + VQ,; = -€i[vi ■ Vvi + V • (Ve, Ve, ; )], 
V • v t = 0, 
d t e t - Aei = Ei[\Vei\ 2 di - Vi ■ VeJ. 

It follows from (|4~T3|) and (|4J5]) that for any 9 € (0, §), 

I 1 - 1 

> Vl|3 ) 3 + (/> l0i|l ) 3 + (X ' Vei1 
and 

/, i "' l f + (^/, ia|i ) l+ (^/,H J ^- 

Applying the WV -estimate to the equation (|4.16[h . we have that V 2 e^ € L 2 (Pi) and 



v2ei ILi(p 7) ~ (iHlWo + HVei||!3 (Pl) ) < c. 



ft 

By the Fubini Theorem and (|4.19[) . we may assume that 



f \V 2 e^<c[ \W 2 ei \i<C. 

JdB,x\-($) 2 ,0] JP 7 



/as3x[-(}) 3 ,o] ip z 
Since (ui, Qi,di) satisfies (|4.6p in Pi, by choosing suitable test functions (j) we have that 

sup f (\ Ul \ 2 + \Vd l \ 2 )+ f {\V Ul \ 2 + \Ad l + \Vd l \ 2 d l \ 2 ) 

-(l) 2 <t<° jB i Jp i 

<C j (\ Ul \ 2 + |VcM 2 ) + (\Qi\ + \ Ul \ 2 + \Vdi\ 2 )\ui\. 
J Pi 

Rescahng (|4.21[) . applying (|4.17|l . and using the Holder inequality, we have 

sup / (H 2 + |Ve 4 | 2 ) + / (|Vw,| 2 + |Ae, -l-ejlVeil 2 ^! 2 ) 



< C / (H 2 + |V ei | 2 ) + (\Qi\ + e ,H 2 + e.lVe,! 2 )!^! < C. 

Jpi 

By the M /2 ' 2 -estimate 14. 8|) and the Sobolev inequality, we have 

/ B jV 2 ei | 2 <P 2 (3)/ B JAe l | 2 + C||Ve l ||^, 2(aB3) 
II 5 

<2? 2 (3) f lAe.p + CHV^II 2 3 

Jb 3 w ^( dB i) 

~& 4 

so that, by integrating over t G [— (-) 2 , 0] and applying (I4.20[) . it holds that 



f |V 2 ei | 2 < 2? 2 (3) [ \A ei \ 2 + c[° ||V 2 ei || 2 

Jps Jps J-ar v 

4 

< C + 2? 2 (3) f |Ae' 2 

JPa 



ll> W 2 -?(dB 3 ) 
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By the point-wise identity | Ae^| 2 = |Ae, + ei\Vei\ 2 di\ 2 + ef| Ve^ | , we have 

f |Ae 4 | 2 =/ |Ae*+e i |Ve i | 2 d*| a + ^ /* |Ve,| 4 . 
J Pi jp 3 J p 3 

3 3 3 

By the Holder inequality, the Young inequality and the Sobolev inequality, we have 



HVe.ll 

L 4 (B 3 ) — II vc «llL 3 (B 3 )H vc 'llL 6 (S 3 ) 



< llVe,- 



lVe, ; 



< 



< 



l|Ve i ||| 3(jB3) (||Ve i - (V e . i )a||| 6(B3) + IKVe^jH 

3 V ~R 



(l + 5o) 2 C 2 (3)||V et ||i 3(B3) ||V 2 e l ||i 2(B3) +C( ( 5o)||Ve l ||i3 (B3) ||Ve l ||i 2(i3 



L 6 (B± 



2 ^WVetWl 

3 ' f 



where (Ve,)3 is the average of Ve. over Bi . Integrating (I4.25|) over t € [— (-) 2 , 0] yields 

4 4 A 



|V ei | 4 < (1 + ^(3)1^*11^^,)/ |V 2 e ,| 2 



+C(<J„)( sup / |VdiH||Vei||£, (Ps) 



-(ir<t<aJBs 



< C(<5o) + (l + ^C 2 (3)||V*||i rL3( p 3) /_ |V^e^ 
Inserting the estimate ()4.26|) first into (|4.24|l and then (|4.23j) . we obtain 



3 Jp 3 

3 



l-(l + <5 ) 2 C 2 (3)2? 2 (3)||Vd l || 2 rL 3 (Pl) / |V 2 ei | 2 



p 3 

3 



< C(<S ) + C / |Ae 4 + e, ; |V ei | 2 d 4 | 2 < C{S Q ). 

JPs 



Therefore, by applying (g^H to we have 



|V 2 ei | 2 < C(5 ). 



Combining the estimates (|4.22j) and (|4.28[) . we obtain 



sup 



(M 2 + |v ei | 2 )+ / (|v Wl | 2 + |v 2 ei | 2 )<c. 



2*2 2 

We may assume, after taking possible subsequences, that 

Qi -> Q weakly in L^(Pi), 

— >■ v strongly in L 2 (Pi), Vv; — > V« weakly in L 2 (Pi), 
ej -> e and Ve, -> Ve strongly in L 2 (Pi), V 2 ei V 2 e weakly in L 2 (Pi). 

Sending i to oo in the equation (|4. 16|) yields that (v, Q, e) satisfies in Pi 

d t v - Av + VQ = 0, 
V • v = 0, 
d t e - Ae = 0. 

Using the Sobolev inequality and interpolations, we see that (|4.29l) gives 

I M 3 + |Q|* + |v e | 3 <a. 
Jpi 
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Hence, by the standard estimates on the linear Stokes equation and the heat equation, we have that for any 
6 G (0, ±), it holds 



(M 3 + |Ve| 3 ) < CO 3 / {\v\ 3 + |Ve| 3 ) < C9 3 



0~ 2 f \Q\^<Ce[ \Q\i <C6. (4.32) 
Jp s J Pi 



In order to reach a contradiction, we need to show that (vi, Qi, ef) converges to (v, Q, e) strongly in L 3 (P2 ). 
To do so, we recall the following Lemma (see [26] ). 

Lemma 4.6. Let X C X c X\ be Banach spaces such that X is compactly embedded in X , X is continu- 
ously embedded in Xi, and X , X\ are reflexive. Then for 1 < a , a\ < oo, 

|u G L a "(0,T;X ) : d t u € L ai (0,T; X x )} is compactly embedded in L a °(0,T;X). 

Now we have the following claims. 
Claim 1. Vi — > v strongly in L 2 (P|). From (|4.29[) and interpolation inequalities, we have 

Ikl^H! + ||Ut||ifiJ(Pi) + ll V ^lli 2 (Pl) < C > 

3 2 2 

1^11^ + \\Ve t \\ L?Ll[ p L] + ||V 2 e t || L 2 (Pi) < C. 

3 1 2 2 

So by the Holder inequality, we have 

k • VVi|3 I |r,|- 



Pi 
7 







)' 













V 3 

These inequalities imply 



|v-(v e ,0v ei )|3 < |v 2 ei | 2 n 



I Veil* <C. 



<H [vi ■ Vvi + V • (Ve* Ve*)] 



"(Pi) 

2 

By (|4.33p and the W^^-estimate of the linear Stokes equation, we have 



< C. 



<9 t u. 



< C. 



i^(p 2 ) 



(4.33) 



(4.34) 



Hence {wj} is bounded in 

Xi = {ue L 2 t Hl{Pz) : d t ueL}Li(P % )\. 
Since Xi is compactly embedded in L\L 2 x (P?f) by Lemma FTol we conclude that Vi — > v strongly in L 2 (P2). 
Claim 2. Ve^ — > Ve strongly in L 2 (P2_). Using (|4]) and the Holder inequality we have 



so that 



Iwi-VeJI 20, < \\vi\\ w, l|Vej|| r4 /p \<C, 



|Ve,;| 2 (ii + Vi ■ VeJI . 30 , < C. 



'iTTCPl) 
3 



Hence the W^-estimate for the heat equation implies 



IftVejll 20 _, 20 < C. 

1 



(4.35) 
(4.36) 



By ([4]) and (|4.36jl . we have {Ve^} is bounded in 

( 20 _ -1 20 ^ 

X 2 = {u G LlHliP*) : ftu G V W a ' (P|)| , 

and so by Lemma [4*^1 we have that Ve^ — > Ve strongly in L?(P£). It is easy to see that by interpolations, 



the claims imply that 



Vi — > v, Ve^ — » Ve strongly in L 3 (Pz 



(4.37) 
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From (|4.37p and (|4.32[) . we conclude that for any 9 G (0, \) and i sufficiently large, 



21 kl 3 + |V ei | 3 <6- 2 I \v\ 3 + |Ve| 3 + o(l) < C6 3 . 



'Pa JPv 

Finally using the estimate (|4.40l) below, with r = 9 and r = i, we have that for any < 9 < I, 



(4.38) 



2 ' |P<|* < c 

'Pe 



\ui\ a + \vdin + e \Pi\i 



After scaling, this implies that for any < 6 < \. 



e- 2 / \Qi\~* < c 



V / (kl 3 + |ve, 
'pi 



Pi 



\Qi 



<C(efe~ 2 + 9). 



(4.39) 



Combining (|4.38p and (|4.39p . we have that for sufficiently large i — i(9), 



(H 3 + |V ei | 3 + |Q 4 |2) <C9. 



Pg 



This contradicts (|4.18[) . if we choose 9 G (0, 4) sufficiently small. 

The next Lemma gives the estimate of pressure function, which is needed in the proof of Lemma 14.5 



□ 



Lemma 4.7. Suppose that (u, P, d) is a suitable weak solution of on P±. Then for any < r < 1 and 

t G (0, j), it holds that 



T JPr 



— r / |P|5 < C 



r\ 2 1 



Pr 



x rW | 3 + |Vd| 3 ) + (^)l^|P|i 



(4.40) 



IB. 



where u r (t) = / u(x,t) for —r 2 <t<0. In particular, it holds that 



r\ J b 



' ' 'P|i<cQ 2 [ sup i 



JP; 



|Vd : 



C 



-r 2 <t<0 r JB r J \ r JP r 

r\ 2 If, .,„ /r\ 1 



(4.41) 



Proof. By scaling, it suffices to consider the case r = 1. Using the equation ([Lip ?, we have 
divdiv [(u - ui(i)) (g) (u - u x (t))] = VjVi ((u - ui(i)) l (u - ui(t)) J ') 

= V, ((« - ui(t)) l V 4 (w - «i(t)y) = V, ((u - UiWjVitt 7 ') 
= Vj(u - ui(f)) l V^' + (u- ui(i))*V i V J -M j 
= (VjU*)(Vitt J ') = VjV i (u i u j ) = divdiv(u0 u). 
Taking the divergence of (|l.l[h . this yields 

AP = - divdiv [(u-ui(i)) ® (u-Mi(i)) + Vd© Vd] . (4.42) 

Let T] G C^°(R 3 ) be a cut-off function of Si, i.e. < 77 < 1, r? = 1 on Si , 77 = outside Si, and |Vr/| < C. 
Define P by 



P{x,t)=- / V^G(x-y) :?7 2 ( 2 /)(( M -ui(i))<»(M-wi(i))+Vd0Vd)(y,i), 
where G is the fundamental solution of the Laplace equation on R 3 . We have 



AP = div div ((u - u x (t)) <8> (u - Ui(t)) + Vd Vd) 
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By the Calderon-Zygmund L p -theory we have 



if) 



rai 5 < / m 



I 5 < 



?? 3 |(tt - (g) (u - ui(t)) + Vd0 Vd| ! 



< 



Si 



(|«-tti(t)| 3 + |Vdl s 



Integrating this inequality over t £ (-r 2 ,0) yields 

C 



\P\2 < 



(| U - Ml (t)| 3 + |Vd| E 



(4.43) 



Pi 



Since the function Q := P — P £ L2 (Pi) satisfies 

AQ(t) = i 

we have by the Harnack inequality that for any < t < i, 



AQ(i) = in Pi, V i G [— ,0] 



3/ 101* <Cr 

T -/Br 



< C7 



IPI 



Integrating this inequality over t G [— t 2 , 0] implies 



< Cr 






f \p\^ + 


f \p\^ 




J Bi 


'Bi J 


+ / 


\u- Ul {t)\ 3 


+ |W| 3 


JB! 








u- Ul (t)\ 3 - 


f |Vd| 3 









(4.44) 



It is now readily seen that (|4.40l) follows by adding the inequalities (|4.43|) and (|4.44p . Using interpolation 
and the Sobolev inequality, we have 



\u- Ul \ 3 <C( / |«| 



B 1 



IVd' 



Inserting (jOSJ) into (jOOJ) yields 

Continuing to iterate the above process, we have 



(4.45) 

□ 



Corollary 4.8. Under the same assumptions as Lemma \4-5\ there exists a £ (0, 1) such that for any 
z\ £ Pro (20) and < r < r < ^P-, it holds 



pah) 



Px(*l) 



|Vd| 3 



< 



(- 



I 



1 



IPI 5 



I 



|vd| s 



(4.46) 



Proo/. Set n = ^ and e a = 2s e . Then it follows from (gj} and ([47X0]) that for any zi G Pra^o), both (|4~9 
and (|4.10p also hold for (u,P,d) with zo,ro and eo replaced by zi,ri and ei respectively. For < p < 7*1 
define $(/o) by 



Hp) 



1 



p ^Pp(^i) 



1 



IPI« 



1 



p ^Pp(zi) 



|Vd| J 



Then applying Lemma T4.5I for (u,P,d) on P ri (zi), there exists 6*o G (0, h) such that for any < r < r\, it 
holds that 

<f(M < i$(r) < 

Iterating ((4]) fc-times, fc > 1, yields 

$(0§r) < 2~ fc $(r). 
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It is well known that this implies that there exists a £ (0, 1) such that for any < r < r < n, $(t) < 
(-) Q $(r). Therefore (TOo) holds. □ 



Proof of Theorem 14.41 We will now prove the smoothness of (u, d) in Pni(zo) by the estimate (I4.46|) . 
The idea is based on the Riesz potential estimates between Morrey spaces, that is analogous to those of 
Huang- Wang [TUJ and Lin- Wang [22] . 

First, let's recall the notion of Morrey spaces on R 3 x R, equipped with the parabolic metric 6: 



s((x,t),(y,s)\ =max{\x-y\,y/\t-s\}, V (x,t), (y,a) 



€ R 3 x 



For any open set U C R 3+1 , 1 < p < +oo, and < A < 5, define the Morrey Space M P > X (U) by 



M*>\V) := \v € Lf oc (U) : \\v\\ p MP , HU) = sup r*~' 



to r < oo 



(4.47) 



zeU,r>o Jp r (z)nu J 
By Corollary 14.81 we have that for some a £ (0, 1), 

u, Vd G M 3 ' 3(1 - Q) (P^.{z j) . (4.48) 

Write the equation as 

d t d -Ad = f, with / := (\\7d\ 2 d - u ■ Vd). (4.49) 

By (|4.48j) . we see that 

/ € Afi' 3 ( 1 - Q ) (Pm(2 )) . 

As in E22 and [10 ], let r) € C 3O (R 3+1 ) be a cut-off function of Pm (z ): < r) < 1, 77 = 1 in Pm(zo), and 
Iftj/I + |V 2 ?7| < CY^ 2 . Set u> = 77 d. Then we have 

d t w -Aw = F, F := ?7 2 / + (d^ 2 - A?7 2 )d - 2Vr/ 2 • Vd. (4.50) 

It is easy to check that F € Mi' 3 ( 1_tt '(l 3+1 ) and satisfies the estimate 



F 



3 



< c 



l + ll/ll M |.3 (1 - Q , (Piii(zo)) 



<C(l + e ) 



(4.51) 



Let T(x,t) denote the fundamental solution of the heat equation on R 3 . Then by the Duhamcl formula for 
(|4.50l) and the estimate (see [TU] Lemma 3.1): 

ivrifot) < x4 „ I fn nNN , v(x,t) ? (0,0), 



<5 4 ((M),(0,0)) : 



we have 



Vw(x,t)\< / |Vr(a:-y,i-s)||J , (tf ) «)|<C 



Jk 3 



S\(x,t),(y,s)) 



where Ip is the Riesz potential of order /3 on R 4 (/3 € [0, 5]), defined by 



Ma) 



\9{y,s)\ 



5{{x,t),{y,s)f-^ 



.9 G L p (R 4 )- 



^C/idJ'D^t), (4.52) 



(4.53) 



3(1 — Q] 0/1 _ \ . 

Applying the Riesz potential estimates (see [10] Theorem 3.1), we conclude that Vic € M !- 2 ° ' 1 QJ (R ) 
and 



V; 



< 



M 1-2" 



< 



Af3' 



l + ll/ll M f.3 (1 - Q ) (Pi!i(2o)) 



<(l+eo). 



(4.54) 



Choosing a t 5 and using lim a ^i j_ 2 ^ = +00, we can conclude that for any 1 < q < 00, Vm € L q (P ro (zq)) 
and 



Vw 



< C(q,r ,e ). 



Since (d — iu) solves 



Li(P r() (z )) 

d t {d - w) - A(d - w) =0inPm(^o), 



(4.55) 
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it follows from the standard estimate on the heat equation that for any 1 < q < +00, Vrf G L q {Pzo (zq)) and 



< C(q,r ,e ). 



i9(Pro (20)) 
4 

Now we proceed with the estimation of tt. Let v : R 3 x [0, +00) — > R 3 solve the Stokes equation: 

' d t v - Av + VQ = -V ■ [r? 2 (V<i 0V(J+u«m)] in R 3 x (0, 00), 
V-w = in R 3 x (0,oo), 

u(-,0)=0 inR 3 . 

By using the Oseen kernel (see Leray [II]), an estimate for v, similar to (14.521) , can be given by 

Hx,t)\<c f f ^r%+i <cii(\x\)(x,t), (MleR^^H, 

Jo Jrs o{{x,t), (y,s)) 3+1 
where X = n 2 (\7d \7d + u ® it). As above, we can check that X G M f^ 1 "") (R 4 ) and 



< C 



ll^llM 3 . 3 ( 1 - a )(P m (zo)) + H U llM3'3( 1 -°)(Pr tt (^o)) 



Hence, by [TO] Theorem 3.1, we have that v € MT^' 3(1 " a) (R 4 ), and 



3(1-°) 3(1-0.) 



< c 



< c 



||Vd||^3,3(l- Q ) (Pla(zo)) + ||«llftf3.3(l-«)(P m (2a)) 



(4.56) 



(4.57) 



(4.58) 



(4.59) 



Af 3' 3 <'- a) (R<) 

By sending a f 5, (|4.59l) implies that for any 1 < q < +00, v e L q (P ro (z )) and 

r , , <C(g,r ,e ). (4.60) 

Note that (it — v) satisfies the linear homogeneous Stokes equation in Phl(zq): 

d t (u-v) - A(tt-u) + V(P- Q) = 0, V • (u-v) = in Pm(z )- 

It is well-known that (u — v) € L°°(P<-o (zq)). Therefore we conclude that for any 1 < q < +00, it 6 
L«(Pm (z )), and 



L^PrO ( Z0 )) 



< C(q, r ,e ). 



(4.61) 



It is now standard that by (|4.56l) and (I4.6ip . and estimates for the linear parabolic equation and the linear 
Stokes equation, (u,d) <E C°°(Pio (z ), R 3 x S 2 ) and the estimate (|4.11j) holds. □ 



5. Existence of L 3 1oc -solutions and Proofs of Theorem 11.21 
In this section, we will prove our main result - Theorem 11.21 



Proof of Theorem 11.21 First, observe that by the scaling invariance of (jl.l[) , (u,P,d) 



[0,7^ 



S 2 solves (jl.ll) under the initial condition (uo,do) if and only if for any A > 0, (u , P , d ) 



[0,T) 



I x S 2 solves (ll.ip under the initial condition (u^dy). Here 
T x = A -2 T; (u^x^^x)) = {\u {\x),do(\x)) for iel 



and 



(u x (x,t),P x (x,t),d x (x,t)) = (Xu(Xx, X 2 t),X 2 P(Xx, X 2 t),d(Xx, X 2 t)) for (x,t) el 3 x [0,T A ). 
Therefore it suffices to prove Theorem 14.41 for R= 1. We divide the proof into six steps. 
Step 1. Approximation of (uq, do) by smooth initial data. We summarize this step into the following lemma. 
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Lemma 5.1. For a sufficiently small eo > 0, let (uo,do) : M. 3 — > M 3 x S 2 , with uo S iui oc (^ 3 ) divergence 
free and (do — eo) G L 3 (K 3 ) for some eo <E S 2 , satisfy 

\\\(u ,Wd )\\\ Llm <e . (5.1) 

Then there exist a large constant Cq > and 

3 

{(4,4)} c C*°°(]R 3 ,M 3 x S* 2 )n P| (L p (R 3 , M 3 ) x W^'^R 3 , S* 2 )) 

p=2 

swc/i i/iai i/ie following properties hold: 

(i) V ■ M§ = m R 3 /or a// fc > 1. 

(ii) As k — > oo, 

(4,4) -> K, do) and Vd k -> Vd in Lf oc (M 3 ) for p = 2, 3. (5.2) 
(mj There exists fc > 1 suc/i £/ia£ /or any > fco, 

|||(4,V4)||Ua (M3) <C eo- (5.3) 

We assume Lemma UTTl for the moment and continue the proof of Theorem ll.2l By modifying the proof of 
the local existence Theorem 3.1 of Lin-Lin- Wang we can conclude that there exist < Tk < +oo and 

smooth solutions (u k , P k ,d k ) : R 3 x [0, r t ]->8 3 xlx S 2 of (fL~T|) . under the initial condition (u k , d k )\ t=0 = 
(4, do). Observe that by applying the proof of Lemma l4~2l with 0=1, the following energy inequality holds: 

{\u k (t)\ 2 + \Vd k {t)\ 2 ) + 2 f [ \Vu k \ 2 + \Ad k + \Vd k \ 2 d k \ 2 = f (|4| 2 + |V4| 2 ), < t < T k . (5.4) 



In particular, we have that (u k ,d k ) £ C([0, T k ], L 2 (R 3 ) x W 1 - 2 { 
Step 2. Uniform lower bounds ofT k . To see this, we first need to show 

Claim. There exists tq > such that if T k is the maximal time interval for the smooth solutions (u k ,d k ) 
obtained in step 1, then T k > to, and 

sup \\\(u k (t)yd k (t))\\\ 3 {m <2C 3 e 3 . (5.5) 
o<t<T i v 

To see (|5.5|) , note that (|5.3|) implies that there exists a maximal time t* k € (0, 7\] such that 

sup |||(w fc (t),Vd fc (t))||||3 (K3) <2C 3 e 3 . (5.6) 

Hence 

|||(^(4),Vd fe (4))lll 3 L3( R 3 ) =2C 3 e 3 . (5.7) 

2 

By a simple covering argument, we see that (|5.6p implies 

sup sup/ (|w & (t)| 3 + |VcZ & (^)| 3 ) < Ceg. (5.8) 

For any fixed x a e K 3 , let O € C^QR 3 ) be a cut-off function of Bi(xq): 

< 0o < 1, 0o = 1 on Si(^o), 0o = outside Bi(xo), and |V0o| < 4. 
For convenience, we set for < t < t* k , 

£ k {<j> 0] {x Q ,t)):= [ [\u k (t)\ 3 + \Vd k (t)\ 3 ]<f> 2 . (5.9) 



2 For K > and < a < 1, first choose the solution space 



X T = {(«,il):l 3 x[0,T]-»l 3 xI 3 :V'ii = 0, V 2 f,dtf 6 C b (R :i x [0,T]) n C a (R 3 x [0,T]), 

(it,d)|t=o = («o>4). IK lt -' u o> £i - <i o)[lcrS' 1 (R8x[o,il) - K }' 
then follow the fixed point argument as in I21| with slight modifications, one can obtain the local existence of smooth solutions. 
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Then by (|3l6]l and (|5~8|) we have that for any < t < t* k , 
d_ 
dt 



-S^ Q] {x ,t)) + {l-Cel) I [\W(\u k (t)\^ )\ 2 + \W(\Wd k (t)\^ )\' 



(5.10) 



< C / (\u k (t)\ 3 + |Vrf fc (i)| 3 )|V0 o | 2 + C sup / \u k {t)\ 3 + \Vd k (t)\ 3 

J»3 yeR 3 \J Bl (y) J 

< Ce 3 + Ce b < Cel. 
Integrating f|5 . 10[) with respect to t £ [0, t* k ) yields 

£ k ^{^tl)) + ^-Cel) f" f [|V(| W fc |i0o)| 2 + |V(|Vd fe |i0 o )| 2 ] 

Jo Js.3 (5.11) 

< Ce 3 t* k + £ k (fo; (a*, 0)) < Ce 3 ^ + C 3 e 3 , 

where we have used (|5.3[) in the last step. Therefore if eo > is chosen such that 1 — Ct\ > 0, then (|5.11[) 
implies 

£|(0o;(xo,4))<^ + C o 3 e 3. 

Taking the supremum of £ k ((f>o; (xo,t k )) over io £ K 3 , we obtain 

2C 3 e 3 = \\\(u k (t* k ),Vd k (tl)\\\ 3 {m < sup E k (<f> ;(x a ,ti)) < Ce 3 t* k + C 3 e 3 . 

This clearly implies that there exists tq > such that Tfc > t* k > To. By the definition of t* k , we also see that 
the estimate (15.51) holds. 



Step 3. Uniform estimation of (u k ,d k ). Note that P k satisfies 

AP k = -div 2 (u fc (g)u fe + Vd fc Vd k ) in R 3 . 
It follows from (|5.4|) . (15"3)) and Lemma I5T2"! that 

sup sup ||P fe (t)-c^(t)|| L3 <Ce , (5.12) 

where c^(t) € R depends on both x £ R 3 and t £ [0,r ]. By (15.51) and (I5.12|) . we see that for any x £ R 3 , 
(u k ,P k — c k a ,d k ) satisfies the conditions of Theorem 14.41 in P^(x ,t ) :— B^(xo) X [0,r ]. Hence by 
Theorem SU we obtain that (u k ,d k ) £ C°°(R 3 x (0, r ),R 3 x S 2 ), and 

su P ||( U fc ,Vd fc )|| cm(R 3 x[5To]) < C(m,S,e ) (5.13) 
holds for any < 5 < ^ and m > 0. 

Step 4. Passage to the limit. Based on the estimates of (u k , d k ), we may assume, after taking subsequences, 
that (u,d) £ P| C^°(M 3 x [<5,t ],R 3 x S 2 ), with (u,Vd) £ L°°([0, r ], L 3 loc (R 3 ), such that 

0<<5<t 

{u k , Vd k ) -> (tt,Vd) weakly ini 3 (R 3 x [0, r ]), (w fc , d fc ) -> in C m (B R x [5,t ]), V m > 0, i? > 0, (5 < r . 

Sending fc — > oo in (|5.8[) yields 

sup || (u,Vd) \\ L 3r m < Ce . 

0<t<T 1 

We can check from (jl.ip and (|5 ,8[) that for any R > 0, 

IK^'^' £ )ILl ([ 0,. o] ,V,-l(B,))^^<+ 00 - 

This implies that 

Vd(i)) ->■ (u , Vd ) strongly in i 3 oc (R 3 ) ast|0. (5.14) 
In particular, we have that (u , Vd ) G C° ([0, r ], L 3 loc (R 3 )). 

Step 5. Characterization of the maximal time interval To. Let To > to be the maximal time interval in 
which the solution (u, d) constructed in step 4 exists. Suppose that To < +oo and (|1.8[) were false. Then 
there exists ro > so that 

limsup|||(u(i),VeZ(t))||| L 3 (R 3) < e . 
t\T 
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In particular, there exists n G (0, ro] such that 

sup \\\(u(t),Vd{t))\\\ L 3 (R 3) < e - 

T -rj<t<T 

Hence by Theoremgai we conclude that (u, d) € C fc °°(R 3 x [0, T ]) n L°°([0, T ], L 3 loc (R 3 )). This contradicts 
the maximality of Tq. Hence (|1.8j) holds. 

Step 6. Uniqueness. Let (iti,di), (u 2 ,do) : K 3 x [0,Tq] -} I 3 x S 2 be two solutions of (jl.lj) . under the same 
initial condition (iio,do), that satisfy the properties of Theorem 11.21 We first show (m,di) = (1*2,^2) in 
R 3 x [0,To]. This can be done by the argument of [27] page 15-16. For convenience, we sketch it here. 
Set u = Ui — u 2 , d = d\ — d 2 . Then (11, d) satisfies 

d t u - Ait = -PV • [m ui - u 2 u 2 + Vdi Vdi - Vd 2 Vd 2 ] 
d t d- Ad = -(m • Vdi - it 2 • Vd 2 ) + |Vdi| 2 di - |Vd 2 | 2 d 2 
(u,d) \ t=0 = (0,0). 

By the Duhamcl formula, we have 

J u(f) = -V[ui iti - u 2 u 2 + Vdi Vdi - Vd 2 Vd 2 ] 
I d[t) = -S[(ui • Vdi - u 2 • Vd 2 ) - (|Vrfi| 2 di - |Vd 2 | 2 d 2 )], 

where 

S/(t) = f e-^'^fis) ds, V/(t) = f e~ (t ~ s)A PV • /(*) da, V/ : R 3 x [0, +00) -> R 3 . 



Recall the three function spaces used in [57J. Let X To denote the space of functions / on R 3 x [0,tq] such 
that 

lll/IHx T0 := sup ||/(t)|U» (ffi3) + ||/||x. < +oo, 

0<t<T O 

where 

||/||x T0 := sup Vt||V/(i)||L»(R3)+ sup (r- 3 f |V/| 2 )^, 

0<t<r xeR 3 ,0<r< v ^i" JP r (x,r 2 ) 

Y T0 denote the space of functions g on R 3 x [0, tq] such that 



IMIy T0 := sup %(t)IU°°(R3) + sup r / |g| < +00, 

0<*<t xeR 3 ,0<r< v /7^ J P r (x,r 2 ) 

and Z ro the space of functions /i on R 3 x [0, ro] such that 

\\h\\z ro := sup VtWhit^L^^ + sup (r^ 3 / |/t| 2 )5 < +00. 

0<t<r o xeR 3 ,0<r< v /r^ JP r {x,r 2 ) 

Since (ti^d*) G i°°([0, r ], L 2 (R 3 ) x M^ 2 (R 3 )) satisfies JET]) for i = 1,2, Theorem 03] and the Holder 
inequality imply that u, € Z To , dj € X ro for i = 1,2, and 

2 

]T(IMIz T „+l|d*llxJ<Ceo. 
i=i 

It follows from Lemma 3.1 and Lemma 4.1 of [27] that 



Mk + ||M|||x T0 < 



(|"i| + |" 2 |)H + (|vdi| + |vd 2 |)|vd| 

+|||u||Vda| + |«i||Vd| + (|Vdi| + |Vd 2 |)|Vd| + |Vd 2 | 2 |d| 

< LE(IKI| X . + lklk )]|Mk + [£(lklk + ll^llx T0 

< e [||u|| ZTO +|||d|||x To ]. 

This clearly implies that (ui, di) = (u 2 , d 2 ) in R 3 x [0, ro]. Since (ui, di) and (ii 2 , d 2 ) are classical solutions of 
(jl.ip in R 3 x [ro, To), and (iti, di) = (u 2 , d 2 ) at t = ro, it is well-known that (u\,d\) = (112, d 2 ) in R 3 x [ro, To)- 
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The proof is complete. 



□ 



Finally, we provide the proof of Lemma 15.11 
Proof of Lemma [SUJ Let 6 G C°°([0, +oo)) be such that 

6(r) = 1 for < r < 1; < 0(r) < 1 for 1 < r < 2; 0(r) = for r > 2. 

Let ?7 g C^°(R 3 ) be a standard mollifier, and define for k > 1 

?yi (x) = k 3 i](kx) and 6> fe (a;) = 6* ( — ) for x £ R 3 . 
k k 

Step 1. Approximation of do- This will be done by two rounds of approximation. It follows from (do — eo) G 
L 3 (R 3 ) that there exists ko > 1 such that for any k > k , it holds 



/ 

Jr 3 \b 



\do -e | 3 < 4- 

By the Fubini theorem, we may assume that for k > fc , it also holds 

'f 9 B k Mo - e | 3 dH 2 < 2/ R 3 VBfc i Mo - e | 3 < 2e 3 , 
' sn PxedB k JdB k nB 2 (x) \^do\ 3 dH 2 < 4 ||Vdo|| L 3( R 3) 
J dBk |Vd | 3 dH 2 < 2j Bk+i \Wd \ 3 < k 3 ||Vd || 3 L 3 (R3) < k 3 e 3 . 

Define the approximate sequence d§ : R 3 — > R 3 by 



d (x) 



if Id < k 



d*(ar) = { (\x\ - k)e + (*+!- \x\)d Q (k-fa) H k < \x\ < k + I 



. eo 



if \x\>k + 1. 



Then by direct calculations we have that 



W fc ll p 
I vu oIIl(b3) 



ivd r 



iv4r 



< 



/ |Vd r+ / |Vd rd^ 2 + / Ido-eol^dif 

7s fc JdBk JdB k 



< k p el < +oo, for p = 2, 3, 



Vdg 



L?(R 3 ) 



~ ll Vd 0|| L 3 (R 3) 



+ sup 

xeaSfc JdB k nB!(x) 

< Ce 3 , 



sup 

x£dB k JdB k C\B 1 (x) 

\Wd \ 3 dH 2 



\d Q ~e \ 3 dH 2 



and for any x G B fc+:L \ B fe 



dist(^(x ),5 2 ) < 



1| JB 1 (x ) 



d o( x o) - <k{v) 



< 



< 



Bi(x ) 



(|x | - fc)e + (fc + 1 - |afo|)do(A!i — r) - do(y) 



\%o\' 



\do(y) - e | 



Bi(io) 



do(y) - do(k-^-) 
\x \ 



(5.15) 



(5.16) 



This implies 



< 



Mo-eo| 3 +||Vdo|| L 3 (R3) <2e . 



sup dist(dg(a;o),5 ) = sup dist(d§(a;o), S 2 ) < 2e 

x eR 3 x <£B k + 1 \B k 
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so that if eo > is chosen sufficiently small then d§(x) remains close to S 2 uniformly for idl 3 . Therefore 



we can project onto S 2 to get cIq(x) 



\d k (x)\ 



for x G R 3 . It is easy to see that dk 



S 2 satisfies: 



d k = d in B fc , d k = e in R 3 \5 fe+ i, 



Vd k 



< Ceo, and 



v4 



For any Z,fc > 1, define d£ ,( (x) = (r?i * dgj (x) for a; G R 3 . Then d G C 



Vd 



k,l 



< Ceo, and 



Vd; 



k.l 



and by the modified Poincare inequality it holds that 



sup dist(do' ; (x),S 2 ) < 

xem. 3 



vd; 



k.l 



< Ce , V I > 1, 



< CFeg < +oo (p = 2,3). (5.17) 
satisfies 

(5.18) 
(5.19) 



< Ck p el < +oo, V I > 1, (p = 2, 3), 



and for any fc > 1, 



lim 



(ll4' i - d o||L ! '(i3 fc _ 1 ) + ll V ( d o''- d o)||L(ij fc _ 1 )) =0, forp = 2,3. 



Therefore, by the Cauchy diagonal process we may conclude that, after taking possible subsequences, there 
exist Ilk) — > oo as k — > oo such that 



d k (x) 



,k,l(k) 



k,l(k) 



(x), V x G R 3 , 



satisfies the desired properties of approximation: dp G C°°(R 3 , S 2 ) n ty 1,p (R 3 , S 2 ) (p = 2, 3), and 

||V4|| L 3 (R3) <C eo, (5.20) 

and for any < R < +oo, 

lim [\\d k o-do\\ LPiBn) + \\V(d k -do)\\ LP(BR) ] = 0, forp = 2,3. (5.21) 

k— >OD 

Next we would like to obtain the desired approximation of uq, whose proof is similar to [2] Theorem 1.4. 
For the completeness, we outline the detail below. 

Step 2. Approximation of uq. Let P : L 2 (R 3 ) — > PL 2 (R 3 ) denote the Leray projection operator. For fc > 1, 
define 

u$(x) =F[6 k uo](x), x G R 3 . 

Since 6 k u Q G L P (R 3 ,R 3 ) and P : LP(R 3 ) -> PLP(R 3 ) is bounded, it follows that V • u k = in R 3 and 
G L p (R 3 ) for p = 2, 3. Now we want to show 



< 



« 



and 
Since 



Lf(R 3 ) 

u§ -> m strongly in Lf oc (R 3 ) for p = 2, 3. 



(5.22) 
(5.23) 



u k {x) = (8 k u )(x) - VA- X V • [9 k u }(x), 
and ||^fcUo|Uf(M3) < ||wo||l ; J(r3), it suffices to show 

|VA- X V ■ [9 k uo] 
Set $ = VA _1 V • [0fctto]. Then we have 



< 



"o 



u k {x) = 9(-)uo(x) - $(x), x G R 3 . 
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It follows from V • uo — that we have 

= VA- X V • [9 k u ](x) = VA-^iVOk) ■ u ]{x) 



K(x-y)V9(l)-u (y) 



where K(x) = c 3 t-^ , c 3 

31^1 

follows. It is easy to see that 



while 



~ [ K(x- y)V9{\) ■ u (y) + \f K ( x ~ 2/)V0(f ) ■ u (y) 

k JB k (x) « k J^\B k (x) k 

x), 

is the kernel of the operator VA _1 . We estimate I and II separately 



(x) 

I{x)+II{x), 

1 



1 

jfe 



V6(-)-u 


< C 


U 




Lf(R 3 ) 





c 



\H(x)\<-l \u (y)\<C 



"o 



so that 

II^IIl3(R3) < C ||M0|Il3(R3) ■ 

Combining these two estimates implies f|5 -22[) . 

For any fixed compact set E CM. 3 and x G E, we write 



C3 



x-y y 



k 7 R 3 Via: - y\ 3 \y\ 
= III k (x)+IV k (u ) 



c 3 / y _„j, 



|y|3 ^ 



Since V#(|) has its support in B^k \ Bk, for k sufficiently large we have that 



x-y , j/ 



and hence it holds 



|^fc(z)| < 



\x-y\ 3 \y\ 3 
C E 



C e 

< -p-, for x € E and y £ B 2 k\ B k , 



fc 4 



E>2k\Bk 



\Mv)\ < 



Ce 



Uq 



as fc — » oo, 



while it is easy to bound /T4(uq) by 



\IV k (u )\ < -3 / |«o(y)| < 

K JB 2k 



u 



Hence we may assume that there exists a constant vector c G K 3 , with |c| < C||ito||.L3( R 3), such that 



lim IV k (u Q ) = c. 

k—tco 



Now we define 



3 x 

u%(x) = v%{x) + -P[0(-)c], x G K 3 . 

£ YZ 



Then we have that u§ G L P (R 3 ) for p = 2, 3, and 



< C 



f(R3) 

It is easy to check that for any x G E, if k — > oo then 



"o 



f )c] . 9 (f )c - VA- V ■ [9 (|)c] = 9 (f) c + „(l) + f £ ^V»(|) ■ c |c. 



Therefore, for any a; G E, if fc — > oo then 



ug^-uoW = (e(~)-l)u (x)-®(x) + ~F[6(~)c 



- l)«o(i) - //J fc (aj) - 7F fe (u ) + |p[0(£)c] -> 0. 
k 2 k 
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This clearly implies (|5.23[) . The proof of Lcmma l5.1l is not complete yet, since it§ ^ C°°(1R 3 , R 3 ). To overcome 
this, we mollify Uq to get 

"o''^) = * u k) ( x )> x e«. 3 ,\/ 1 >i. 



Then it is straightforward to check that Uq' 1 € C°°( 



fW(M 3 ,IR 3 ) forp = 2,3, V • u, 



k.i 



0. 



k,l 



L\ (R 3 ) 



< 



L\ (R 3 ) 



< c 



Uq 



and for any k > 1, 

-Up'' — > Uq strongly in Lf oc (R 3 ) for p = 2, 3, as / — » oo. 
Thus, by the Cauchy diagonal process we may assume that there exist l(k) — > oo as k — > oo such that 

4(x) = u k m (x), iel 3 

satisfies the required properties of approximation of uq : Uq € C°° (M 3 , M 3 ) n £ p (R 3 , M 3 ) for p = 2 , 3 , V • Uq = , 

Lf(R3)' 



u k 




< c 


u 


L 3 (R 


3 ) 





and 

Mo u Q strongly in if oc ( 
This completes the proof of Lemma 15.11 



for p = 2, 3, as k — > oo. 



□ 
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